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Abstract 

In the framework of melting as a dislocation-mediated phase transition we de- 
rive an equation for the pressure dependence of the melting temperatures of the 
elements valid up to pressures of order their ambient bulk moduli. Melting curves 
are calculated for Al, Mg, Ni, Pb, the iron group (Fe, Ru, Os), the chromium group 
(Cr, Mo, W), the copper group (Cu, Ag, Au), noble gases (Ne, Ar, Kr, Xe, Rn), 
and six actinides (Am, Cm, Np, Pa, Th, U). These calculated melting curves are 
in good agreement with existing data. We also discuss the apparent equivalence of 
our melting relation and the Lindemann criterion, and the lack of the rigorous proof 
of their equivalence. We show that the would-be mathematical equivalence of both 
formulas must manifest itself in a new relation between the Griineisen constant, 
bulk and shear moduli, and the pressure derivative of the shear modulus. 
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1 Introduction 

The idea that a proliferation of dislocations is associated with melting dates back to Mott 
n. The very first theory of dislocation- mediated melting was a success, inasmuch 
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as it predicted a first-order transition, as a consequence of incorporating the mutual 
screening of dislocations, in agreement with observations. Molecular dynamics 0] and 
Monte Carlo ||J] calculations have more recently provided further evidence for the notion 
that dislocations drive the melting transition in three dimensions. There is also some 
experimental evidence that line defects are present in solids near melting ||. 

In refs. |J [7| we formulated a dislocation theory of melting in which dislocations near 
melt were modeled as non-interacting strings on a lattice. The possible configurations of a 
dislocation were taken to be closed random walks. Screening of long-range strain fields by 
other dislocations in a dense ensemble results in a — p In p dependence of the free energy 
on the dislocation density, p, and thus a first-order transition. We obtained the following 
relation between the melting temperature T m , the shear modulus, G, the Wigner-Seitz 
volume, Vwsi the coordination number, z, and the critical density of dislocations, p(T m ) 
(in units where kg = 1) : 



kXGvws ( <\~ 



Tm 87rln(/- J l) ln U 2 P(T m ))' (1) 



Here b is the length of the shortest perfect-dislocation Burgers vector, k is 1 for a screw 
dislocation and (1 — u)^ 1 « 3/2 for an edge dislocation (v is the Poisson ratio), A = b 3 /vws 
and a, which accounts for non-linear effects in the dislocation core, has a value of 2.9 0. 
Experimental data on 51 elements show that 

GVWS 1.01 ±0.17 (2) 



4irT m \n(z- 1) 

at zero pressure ||. Eqs. (1) and (2) imply that the critical dislocation density at zero 
pressure is 

p{T rn ) = (0.61 ±0.20) ft" 2 . (3) 
This value is in good agreement with the critical density 

p(T m ) = (0.66 ±0.11) r 2 , (4) 

obtained by applying our relation for the latent heat of fusion 0, 

L m = - b 2 p(T m )RT m ln(z - 1), (5) 

to data on latent heats for 75 elements. Hence, b 2 p(T m ) is approximately constant across 
the Periodic Table with the numerical value 

b 2 p{T m ) = 0.64 ±0.14, (6) 

which is the uncertainty- weighted average of Eqs. (3) and (4). 

In this paper we investigate the validity of our melting relation, Eq. (1), up to pressures 
of order 100 GPa, by comparing to experimental melting curves, i.e., melting tempera- 
tures versus pressure, p. This comparison requires vws{p)i or hs equivalent, the pressure 
dependence of the compression, rj = V /V. We obtain rj(p) from the bulk modulus, B{p), 
which is extrapolated to high pressure using only its value and first pressure derivative at 
ambient conditions, viz., room temperature and zero pressure. The shear modulus G(p) 



2 



is similarly extrapolated to high pressure. Pressure derivatives of G and B are typically 
0(1), so the extrapolation of the bulk modulus is expected to break down at pressures 
of order the ambient bulk modulus. The parameter k in Eq. (1), which depends on the 
Poisson ratio, varies by only a few percent between p = and 100 GPa (we discuss this 
in more detail in Section 2). Since the accuracy of our melting relation at zero pressure is 
17%, we take k to be a constant. We also make the necessary but reasonable assumption 
that b 2 p{T m ) is also a pressure-independent constant. With this assumption, we find that 
our melting relation agrees well with experimental melting curves up to pressures ~ B, 
and, in fact, our extrapolation of T m is often in good agreement with data to pressures 
fa 2B. In addition to the good agreement with the existing melting curve data, we also 
predict the high-pressure melting curves of Ag, Au, Cr, Cu, Mo, Os, Ru, W, and several 
actinides. 



2 Melting curve equation 

We now consider the pressure dependences of the factors appearing in our melting relation, 
Eq. (1). The parameter A is constant by its definition, a is also assumed to be a constant, 
and k may be taken as constant provided that the Poisson ratio v has a very weak pressure 
dependence. In fact, for an isotropic medium || 

1 3B - 2G 

v = . (7) 

2 W + G y 1 

Although both G and B vary with pressure, the ratio in Eq. (7) varies only weakly. 
Consider, for example, Cu, for which v fa 0.34 at p = 0. At p = 100 GPa, we calculate the 
values of G and B with the help of Eqs. (13) and (14) below, with their pressure derivatives 
taken from ref. 0, and find v ~ 0.38. Therefore, in this case the corresponding values of 
1/k fa 1 — v/2 [|7J are 0.83 and 0.81, respectively, so that the variation in the average value 
of 1/k is ~ 2%. Thus, the pressure dependence of 1/k can be safely neglected. (There 
exists an upper bound on the change in the value of 1/k with pressure. In the ultra-high- 
pressure limit, p oc r/ 5//3 , in agreement with the theory of the free electron gas (Fermi gas). 
Therefore, B = —Vdp/dV = r]dp/dr] oc rf^ ^> G oc r/ 4//3 , and hence v — ► 1/2, in view of 
Eq. (7) (see also flOl). Thus, in contrast to the Poisson ratio which changes by ~ 50%, 
1/k fa 1 - u/2 changes by fa 10%: 5/6 3/4.) 

We assume further that the mean interdislocation spacing at the melting point, R ~ 
1/ \J p(T m ) , scales with b, independent of pressure, and hence b 2 p(T m ) is a pressure- 
independent constant (with a numerical value of 0.64 ± 0.14, in view of Eq. (6)). It 
then follows from Eq. (1) that, provided the coordination number does not change with 
pressure (i.e., the element either remains in the same crystalline phase, or changes phase 
without changing the coordination number, e.g., a face-centered cubic structure *-> a 
hexagonal close-packed structure), the melting relation is given by 

G(p,T m (p))vws( P ,T m (p)) = congt (g) 

Tm\P) 

The dependence of vws on pressure and temperature is just the equation of state of 
the metal. Let us first focus on its temperature dependence. The fixed-pressure ratio 
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of Wigner-Seitz volumes at T m and T = is equal to 1 + (3T m , where (3 is the volume 
expansivity. At p = 0, j3 is typically of order K~ 4 , and melting temperatures are at 
most about 4000 K, so v ws changes by only a few percent between T = and T m . Assum- 
ing that j3 does not increase appreciably with compression, we can use room-temperature 
values for vws- 

In contrast to VwSi the dependence of G on T is not necessarily weak. Its T-dependence 
involves two characteristic temperatures, namely the Debye temperature, Tjj, and the 
melting temperature. G is always monotonically decreasing with T, and is nonlinear for 
T ~ Td and linear from To to T m . However, there are no experimental data, no computer 
calculations, and no theoretical guidance that tells us how the temperature dependence 
of G varies with pressure. In particular, how does the (negative) slope of the linear 
region vary with pi At this point we have no choice but to conjecture. We assume that 
G(p,T m (p))/G(p,0) is a slowly varying function of p, so it can be considered constant up 
to moderate compressions, say, 20% to 30%. Thus, G(p, T m ) is replaced by G(p, 0) in Eq. 
(8). In addition, data on the p = temperature dependence of shear moduli []TTJ clearly 
show that G(p, 300) ~ G(p, 0), and therefore, we use the room temperature value of the 
shear modulus in our melting relation. 

Subsequently, the explicit dependence of G and v ws on T will be dropped. It will be 
understood that G and vws are at room temperature. Our melting relation now reads 

G( P )v ws (p) = congt (g) 



T m (p) 

Differentiating Eq. (9) with respect to p, one finds 

1 dT m _ 1 dG 1 
T m dp G dp B' 
where we have used the definition of the bulk modulus, 



(10) 



BO,) . -V% = (11) 
dv dvws 



Thus, upon integration, Eq. (10) gives 



t.m_gm r (12) 



T m (0) G(0) r I h S(P' 
To proceed further, we have to specify G{p) and B{p). 



2.1 The shear modulus G at finite pressure 



For the shear modulus at all pressures, we use the relation |L2 

P 

,1/3 



G — G + G' (13) 



where G' = (dG/dp)o. The subscript refers to ambient conditions: T ~ 300 K and 
p = 0. 

This equation satisfies the requirement that G oc 7/ 4//3 as r] —> oo, since p oc i] 5 ^ 3 . With 
the values of G' for 32 elements tested in ref. Jl2[ Eq. (13) gives nearly the right value for 
the proportionality constant between G and r/ 4 / 3 at high compressions. Eq. (13) works 
well for a diverse selection of engineering metals covering many different crystal structures 



and nearly all groups of the Periodic Table |12 
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2.2 Compression and the bulk modulus B at finite pressure 

Expanding the bulk modulus around p = we have 

1 



B(p) = B 



■B p + ..., 



(14) 



where B Q and B' Q = (dB/dp) , Bq = (d 2 B/dp 2 )o, . . . can be extracted from equation of 
state data. Values of Bq are known for a few elements only (their determination is highly 
uncertain and involves an error of order 100% |H|), and besides, B' ' first appears in the 
(p/Bq) 3 term in the power series expansion of rj : 



rj = exp 



1 + 



p_ 

Bo 



p dp' 
o B(p') 



B'-l 



2B + (B' 



B' 2 



2B Bg)p 



2B + (B' 



B' 2 



2B BZ)p. 



l/^/B[?-2B B>> 



P_ 

B 



+ 



BqBq 
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p_ 

Bo 



(15) 



Since only Bo and B' are generally known (for almost all the elements, see ref. 0), 
we restrict ourselves instead to the first two terms in Eq. (14). Then the compression 
simplifies to 

T) = 1+ -^p 



Bn 



(16) 



Eqs. (15) and (16) are two different approximations to the Murnaghan equation of state 
It then follows from Eqs. (12), (13) and (16) that the equation of the melting curve is 



T m (p) = T m (0) M + 



G' 



B' 



-1/3B' Q 



1 + ?rP[ 1 + W p 



(17) 



As discussed in Section 4, this equation is only valid for pressures p ~ 2B. 
It follows from (17) that for p«B 



T m {p) = T m (0) 



BqG'q 

Go 



-'"5 



3 G 



B'n + 1 



For the vast majority of the elements, B' > 5/3 and B' approaches 5/3 in the limit of large 
compressions. (In this limit p oc r/ 5 / 3 , and therefore B = —Vdp/dV = rjdp/dr] = 5p/3, 
i.e., B' = 5/3.) In fact, the average value of B' Q for the 65 elements analyzed in 0, except 
for Ce for which B' < 0, is 4.30 ± 1.40. Hence, if 

G'n 3 B' n + 1 



> 



Go 8 Bq 



(19) 



it follows from B' Q > 5/3 that also G' Q /Go > 1/B , i.e., Eq. (18) is of the form T m (p) = 
T m (0)(l + ap — bp 2 + ...), a, b > 0, and describes melting curves for which melting tem- 
peratures increase with pressure [1(J. If, however, 



Go B 



(20) 
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and B f > 5/3, then also G f /G < 3/8 (B' + l)/B , i.e., Eq. (18) is of the form 
T m {p) = T m (0)(l — ap + bp 2 — •••), a,b > 0, and describes melting curves for which 
melting temperatures initially decrease with pressure [16|. For Si, for example, with the 
data from ref. we find G' /G <1/B and B' = 4.19, in agreement with the negative 
initial slope of the experimental melting curve. Eqs. (19) and (20) plus B' Q > 5/3 should 
be considered our criteria for the two types of melting curves discussed above. 



3 Melting curves: comparison with data 



In this section we compare our melting curve, Eq. (17), to some experimental melting 
curves, and predict a number of melting curves that can be compared with experiment in 
the not- so- distant future. 

We have found 5 elements for which melting curves have been measured to higher 
pressures, p ~ 0(100 GPa): Al, Fe, Ni, Pb and U. We compare experimental data for 
these elements with our curves in Figs. 1-5. For Al, we also show the best fit to data in 
the form of the Simon equation, T m {p) = T m (0)(l + ap) b f2"7|| . For Fe, the experimental 
data are from ref. |28|| , and from ref. |29] for Ni. For Pb, we combine the high-pressure 
data of ref. |$0| with the low-pressure data of ref. |31[ as corrected in ref. |32[]. For U, the 
high-pressure data of ref. |3T 



are combined with the low-pressure data of ref. [34] . 



As claimed in ref. |27|], the Simon equation may not be the best functional form for 
a fit to data. In fact, the initial slope provided by this equation for the Al melting 
curve is 80 K/GPa, in contrast to 59 and 65 K/GPa from the two previous low-pressure 
measurements P7| . This accounts for the difference between the two curves in Fig. 1. 

In Fig. 4, in addition to Fe, we also plot melting curves for Ru and Os, elements in 
the same column of the Periodic Table. Those curves should be considered predictions 
for these metals. 

In Fig. 5, in addition to U for which there are high-pressure data, we also plot melting 
curves for the 5 actinides Am, Cm, Np, Pa and Th. For Np, we also show the low-pressure 
data of ref. |18| . We do not show the low-pressure data of ref. |17] for Am since they would 



overlay the low-pressure U data. We have checked that our melting curve is in agreement 



with the low-pressure Am data. For Cm, the values of Bq and B' are taken from ref. |3l| , 
and the value of Go is that estimated in ref. [Q. For Pa, the values of Bq and Go come 
from ref. ||. We estimate the values of G' for Cm and Pa from Th and U, their neighbors 
in the same row in the Periodic Table. Our earlier G' estimates for Am and Np lead to 
7 = 1.05 and 1.09, respectively, in Eq. (24), which implies that such estimates are reliable. 
The values of B' for Np and Pa are also estimated by interpolating between Am, Cm, 
Th and U. (We note that this estimation of B' Q is justified by the pronounced periodic 

The value of B' Q for Am is taken from |17|| . We emphasize that 



behavior of B' n in Z 36 



the predicted melting curves assume constancy of coordination number along them. In the 
case of Am, e.g., there is still disagreement over the correct sequence of phases and their 
transition pressures ||16[| , so this assumption may well be incorrect. For Th, however, it is 
claimed that there is a transition from a face-centered cubic structure to a body-centered 
tetragonal structure that changes coordination number . This transition occurs in the 



pressure range of 70 — 100 GPa ||37]| , and thus our predictions for the Th melting curve 
up to 75 GPa should be quite reliable. 

Although we can account for a decrease in melting temperature with pressure in our 
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theoretical framework (Eqs. (17), (20)), we do not consider such cases here, among which 
there are Pu and Ce. It has been established that for Th, which is in the same column 
as Ce, AV > 0, and therefore, in view of Eq. (23), its melting temperature increases with 
pressure. 
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Fig. 1. Melting curve for Al. The dashed line is the Simon-fit to the data of ref. 
which are not shown explicitly. The diamonds are the low-pressure data from ref. 
The triangle is the shock-melting point at 125 GPa from ref. [B3. The boxes are the 
points at 25, 69 and 137 GPa calculated in ref. |40| from shock-melting data. They are 
assigned 20% error bars P0|. 
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Fig. 2. Melting curve for Ni. The diamonds (with small error bars) are the data of ref. 
p9| . The boxes are the points at 79 and 250 GPa calculated in ref. |4(| from shock-melting 
data. The corresponding error bars are not quoted in ref. EDI. 
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Fig. 3. Melting curve for Pb. The diamonds are the low-pressure data of ref. [31 corrected 
as in ref. f3"2]| . The triangles are the data from ref. [[HJ, and the dashed line is a best fit 
p0| to the data. The boxes are the points at 12, 34 and 68 GPa calculated in ref. 
from shock- melting data. The corresponding error bars are not quoted in ref. fl40f| . The 



star is the point at 118 GPa calculated in ref. [41] from shock-melting data 
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Fig. 4. Melting curves for the elements of the iron group (Fe, Ru, Os). The diamonds 
are the data of ref. 



28 , and the dashed line is a best fit 



to the data. The boxes are 
the shock-melting points at 235 GPa and 300 GPa El. The triangle is the shock-melting 



point at 240 GPa [43]. The star is the shock-melting point at 243 GPa [44 . 
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In Fig. 6 we plot the low-pressure data of ref. [25] for W, and our melting curves for 
W, Mo and Cr. The initial slope of our melting curve of Mo, 26 K/GPa, is consistent with 
that predicted in ref. ||22|| : (34 ± 6) K/GPa. The same melting curve gives T m ~ 9650 K 
at p = 390 GPa, in good agreement with the shock-melting temperature ~ 10000 K at 



the same pressure, found in ref. [45 




25 50 75 100 

GPa 

Fig. 6. Melting curves for the elements of the chromium group (Cr, Mo, W). The data 
for W are from ref. [EHl. 



In Fig. 7 we compare the low-pressure data of ref. for Cu, Ag and Au with our 
corresponding melting curves. Although the initial slopes of these curves are somewhat 
less than those of the data (the corresponding values of 7 in Fig. 1 are ~ 0.8), they are in 
good agreement with the best extrapolation of data to higher pressures made in ref. [|6| , 
and with the calculation of ref. [^0j in the case of Cu. 
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Fig. 7. Melting curves for the elements of the copper group (Cu, Ag, Au). The diamonds, 
the stars, and the triangles are the low-pressure data of ref. [f5| in the inset, and the 
best extrapolations of these data to 20 and 30 GPa in the main plot for Cu, Ag and Au, 



respectively. The boxes are the points at 45 and 128 GPa calculated in ref. j|0| from 
shock-melting data for Cu. The corresponding error bars are not quoted [HO]. The gray 



diamond is the shock-melting point for Cu at 37 GPa [B7 



In Fig. 8 we compare the low-pressure data on the noble gases to our corresponding 
melting curves. The unknown values of G' for Ne, Ar, Kr and Xe are calculated with 
the help of Eq. (25) below using the measured values of B , Gq and 70 (the Griineisen 
constant) [|7|. The values of B' for Ne, Ar, Kr and Xe are taken from f4~8fl . In the case of 
Rn, for which B , B' Q , Gq and G' Q have not been measured, we first calculate Gq using the 
(approximate) relation GV/T m = const for the noble-gas group, where V = vws^a is the 
molar volume. This relation follows from Eq. (1) provided that n, A, a and z do not vary 
within this group. The value of the constant is determined by using the corresponding Ne, 
Ar, Kr and Xe data in this relation. We then calculate Bo using Eq. (7) with the value 
of the Poisson ratio for Rn determined by extrapolating from the corresponding values 
for Ne, Ar, Kr and Xe. Finally, we determine both B' Q and G f by again extrapolating the 
Ne, Ar, Kr and Xe data. 
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Fig. 8. Melting curves for the noble gases Ne, Ar, Kr, Xe and Rn. The diamonds are the 



data of ref. [49]. The stars are the data of ref. |pOfl . The boxes and triangles come from 



ref. [51 



Finally, in Fig. 9 we show the experimental data and our theoretical melting curve for 
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Fig. 9. Melting curve for Mg. The low-pressure data are from ref. ||52| . The high-pressure 



data are the shock-melting points of ref. [53 
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The 24 melting curves considered above constitute convincing evidence for the validity 
of our formula for melting temperature as a function of pressure, Eq. (17). 



4 The range of validity of the new melting curve 
equation 

In deriving our melting curve, Eq. (17), we have used both Eq. (13) for the pressure 
dependence of the shear modulus and the Murnaghan equation of state, Eq. (16). Since 
Eq. (13) has the correct zero-pressure limit (its Taylor series expansion in p at p = is 
G = G + G' p — (G' /3B )p 2 + . . .) and is claimed to have the correct ultra-high-pressure 



limit | I2| , we assume that this equation is valid over the entire pressure range. In any 
event, we do not have data to either confirm or invalidate this assumption. It then follows 
that the range of validity of Eq. (17) depends crucially on the range of validity of the 
Murnaghan equation of state, Eq. (16). 



The Murnaghan equation of state was examined in ref. |15| , together with a number 
of different equations of state, by comparing with the theoretical results calculated by the 
augmented-plane-wave method and the quantum-mechanical model proposed by Kalitkin 



and Kuz'mina |54| from low to ultra-high pressures. It was shown that the Murnaghan 
equation is in good agreement with the theoretical results up to V/Vq — 0.7, i.e., up to 
compressions ~ 1.4 — 1.5. Since for the vast majority of the elements B' « 5 0, we 
conclude, on the basis of Eq. (16), that the Murnaghan equation, and consequently, our 
equation for melting curve, Eq. (17), is valid up to pressures p ~ 2B . The melting curves 
for Al and Pb in Figs. 1 and 3, respectively, and for Ne and Ar in Fig. 8 show that in 
some cases Eq. (17) is good to pressures even greater than 2B . 

For reliable predictions of melting curves to much higher pressures, p ~ 1 TPa, one 



has to use a better equation of state than Murnaghan's. Hama and Suito [15] claim that 



the Vinet equation of state |55[] is consistent with first-principles theoretical calculations 
to compressions r\ ~ 5. Reference |56| also finds the Vinet equation of state to be most 
accurate among various suggested equations of state. In fact, we have calculated that the 
melting temperatures for Fe and Pb at pressures p ~ 50I?o as given by Eq. (17) are about 
two times higher than those given by a relation that derives from Eqs. (12), (13), and the 
Vinet equation of state. 

Another possible source of disagreement between the new melting curve, Eq. (17), and 
data may be inaccurate values of elastic constants and their pressure derivatives in some 



cases. The Murnaghan equation and its frequently used partner - the Birch equation |h] 
- are derived from the second-order Taylor series expansion of the bulk modulus [as in Eq. 
(14)] or the elastic strain energy with respect to pressure or strain, respectively. Thus their 
validities are, in principle, restricted to a narrow range of compression. Extending this 
range would entail the inclusion of higher-order terms. This could explain why the values 
of B , and especially those of B' and Bq, obtained from experiments which cover different 
ranges of compression by using a fitting method, are usually different. In many cases these 
differences between different experiments are small and can be safely neglected. In some 
cases, however, they are large, and so their use for predicting physical observables, such as 
melting temperature, is dubious. For example, in the case of Ni, we have used the value 
B' Q = 6.20 given in |§. Reference f58|, however, quotes B' ~ 30 (!). Similarly, for Mo we 
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have used B' Q = 4.4 of ref. 0, while ref. |53[] gives B' Q ~ 20. (We note that the use of the 
values B' = 6.20 for Ni and 4.4 for Mo is justified in view of the recent compilations of 
experimental data on B' ||59|| .) Although the numerical value of B' does not matter at 
low p, since it first appears in the (p/B Q ) 2 term, in view of Eq. (15), it would strongly 
affect the predicted melting curve at pressures p ~ O(B ). 

There are also inconsistencies in the values of Go quoted in the literature. For example, 
for Pb we use the value G = 8.6 GPa from ref. |§, whereas ref. |§ quotes G = 5.5 GPa. 
(Our own calculation ||, based on the values of the elastic constants en, c\i and C44, 
shows that 8.6 GPa is preferred over 5.5 GPa.) Likewise the values of Go for K and Na 
from ref. || are 0.9 and 1.98 GPa, whereas ref. || quotes 1.3 and 3.5 GPa, respectively. 



5 Relation of dislocation-based melting relation to 
the Lindemann criterion 

The well-known Lindemann melting rule is based on the assumption that all elemental 
solids melt when the atomic vibrational amplitude is a fixed pressure-independent fraction 
of the interatomic distance. As shown by Lindemann [60], this implies the invariance of 
the Lindemann number 

e D (^-) 1/2 v^ = L (21) 

"111' 

along the melting curve. Here Op is the density-dependent Debye temperature, V is the 
molar volume, and M is the molar mass. It is found that L « 150 ||. 

There are compelling reasons to suppose that our dislocation-based melting relation 
is somehow equivalent to the Lindemann criterion. First of all, Eq. (21) gives melting 
curves that are typically very close to those predicted by our dislocation-based melting 
relation. (For example, the melting curve for Mg in Fig. 6.3 of ref. [0 is very similar 
to our curve in Fig. 9.) Furthermore, the Lindemann number, which is proportional to 
the ratio of atomic vibrational amplitude to the lattice constant at melt, is analogous to 
b 2 p(T m ), since both are presumed constant along the melting curve. In the dislocation- 
based approach, melting is associated with a critical configuration of dislocations, and 
for any such configuration there is a corresponding mean displacement of atoms from 
their equilibrium positions. Hence L and b 2 p(T m ) are clearly related, and therefore the 
left-hand sides of Eqs. (9) and (21) are related as well. 

The mathematical equivalence of our melting relation and the Lindemann criterion 
would be established if it could be determined that the left-hand side of Eq. (21) is a 
fixed fraction of the left-hand side of Eq. (9). A search of the literature has turned up 
two results which show that L 2 is approximately proportional to Gvy/s/T m . For a Debye 
solid the relation is |H] L 2 = f(v(p,T))Gvws/T m , where / is a complicated function of 



v. Thus the two melting relations are not rigorously equivalent. 

A second connection between the two melting formulas is provided by the following 
approximation for the Griineisen constant ||, 
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where 

1S{P) = —^P)-^ (23) 

( ) _ 1 B-( P ) dW + jGjp)) 1 
TlW "2BS( p ) + | G (p) rfp 6 ' ^ 

are the contributions of the shear (transverse) and longitudinal acoustic modes. Here B T 
is the isothermal bulk modulus, which is equivalent to B that we are using in this paper. 
Eqs. (22)- (24) follow from the two assumptions that (i) the only appreciable contribution 
to the heat capacity of a crystal arises from lattice vibrations, and (ii) averaging over all 
modes is equivalent to averaging only over the low- frequency acoustic modes. (I.e., the 
contribution of the optical modes is equal to that of the acoustic modes.) If in addition 
it is assumed that B s (p), the isentropic bulk modulus, is proportional to G(p), then 

/ n G'(p) B(p) 1 , N 

ls(p) = 1l(p) = lip) = -f 1 ^ - g • (25) 

However, there is no basis for this assumption, i.e., jsip) 1l(p) is to be expected. For 
example, in the ultra-high pressure limit, B s (p) ~ B T (p) = 5p/3 and G(p) ~ p 4 / 5 , quite 
different dependencies. 

Integration of Eq. (25), using B{p) = —dp/d\iaV(p) and ^(p) = —d\n9E>(p)/d\iaV(p), 
gives 

»Up)v^p)Itm _ const 

G(p)V(p)/T m (p) C ° nSt ' (A) 

that is, L 2 oc Gvws/T m - We emphasize that this proportionality is founded on a number 
of uncontrolled approximations. 

Equation (25), which would ensure a rigorous mathematical equivalence of the defect 
and mechanical (Lindemann's) approaches to melting, does not follow from first principles. 
This means that the defect and mechanical approaches to melting are basically different. 
Moreover, since the mechanical approach does not have a solid thermodynamic basis, it 
cannot, for example, predict the latent heat of fusion. In contrast, the defect approach 
predicts the latent heat of fusion, Eq. (9), which is in good agreement with data for 
three-quarters of the Periodic Table 0. 

Finally, we wish to make the following comments on Eq. (25). We did not check 
extensively its validity at zero pressure, since that would go beyond the scope of this 
paper. We do, however, have some evidence that Eq. (25) is rather well satisfied: with the 
data from ref. M, we calculate from the above relation 7 = 2.25 vs. measured 2.40 for Ag, 
3.08 vs. 2.99 for Au, 1.66 vs. 1.78 for Fe, 1.28 vs. 1.29 for K, and 1.18 vs. 1.19 for Na. We 
have actually used Eq. (25) in Section 3 to calculate G' Q for noble gases in order to get their 
melting curves via Eq. (17) and to compare with experiment. Good agreement between 
the calculated and experimental curves is another hint on the approximate validity of this 
formula. Also, Eq. (25) has the correct ultra-high-pressure limit in which 7 — > 1/2 [ TCfl . 
since in this limit G ~ p 4 ^ 5 and B = 5p/3. 
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6 Concluding remarks 



We have extended the framework of melting as a string-mediated phase transition to non- 
zero pressure and derived a new equation for the melting curve, Eq. (17). As discussed 
above, with accurate experimental values of all the parameters involved, this equation 
reproduces the existing experimental melting data, and predicts unknown melting curves 
to pressures p ~ 2B . For higher pressures, a better equation of state than Murnaghan's 
should be used, e.g., the Vinet equation of state. 

We have addressed the apparent equivalence of defect and mechanical approaches to 
melting curve, and demonstrated that both approaches are basically different. We have 
shown that their would-be rigorous mathematical equivalence must manifest itself in a 
new relation, Eq. (25), which we have not tested in detail. 

To summarize, we have calculated melting curves for 24 elements: Al, Mg, Ni, Pb, the 
iron group (Fe, Ru, Os), the chromium group (Cr, Mo, W), the copper group (Cu, Ag, 
Au), the noble gases Ne, Ar, Kr, Xe and Rn, and the six actinides Am, Cm, Np, Pa, Th 
and U. These calculated melting curves are in good agreement with existing data. 
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